Thermally- Assisted Current-Driven Domain Wall Motion 
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Starting from the stochastic Landau-Lifschitz-Gilbert equation, we derive Langevin equations 
that describe the nonzero-temperature dynamics of a rigid domain wall. We derive an expression 
for the average drift velocity of the domain wall (raw) as a function of the applied current, and 
find qualitative agreement with recent magnetic semiconductor experiments. Our model implies 
that at any nonzero temperature (raw) initially varies linearly with current, even in the absence of 
non-adiabatic spin torques. 
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Introduction - Possibilities opened up by modern 
nanofabrication capabilities have motivated renewed in- 
terest in current-induced domain wall motion, a phe- 
nomenon first predicted and subsequently observed in 
seminal work by Berger Q, The theoretical El El 
S II II S and experimental [H El H 13 Mm 
study of domain wall motion induced by spin transfer 
torques yj| [XJ3, UM lill is currently one of the most ac- 
tive subfields of spintronics j^. Recent research has 
highlighted a number of fundamentally interesting issues 
that are currently under spirited debate. One contro- 
versy concerns intrinsic pinning 0-01 . i.e., domain walls 
that are stationary up to a critical current in the absence 
of spatial inhomogencity. An intellectually distinct but 
phenomenologically related debate surrounds theories of 
nonadiabatic spin torques, which differ widely in their 
predictions 0, H, IS E| ■ ^ turns out that in the pres- 
ence of these torques, a domain wall is never intrinsically 
pinned 



This Letter is motivated primarily by the recent ex- 
periments of Yamanouchi et al. 15], in which current- 
induced domain wall motion was studied over five or- 
ders of magnitude of average wall velocity. An important 
conclusion of these authors is that at low temperatures 
the domain wall undergoes creep motion, i.e., that the 
domain wall does not move rigidly. Yamanouchi et al. 
arrive at this conclusion because the effects of nonzero 
temperature, when treated |2l) in a rigid domain wall 
approximation, seem to lead to results that are irrecon- 
cilable with experiment. In this Letter we demonstrate 
that a systematic theory of the influence of a thermal 
bath on the current-driven motion of a rigid domain wall 
leads to results that are in qualitative agreement with 
experiment. In the following sections we first explain our 
theory of nonzero-temperature domain wall motion and 
then discuss its implications for recent experiments. 



Drift velocity of a rigid domain wall — Our starting point 
is the stochastic Landau-Lifschitz-Gilbert (LLG) equa- 
tion [H HI |H |H HI IHIH for the direction of mag- 



netization Q 



(1) 



where H is the effective field defined by H(x) = 
— <5£?MM[^]/(Mf2(x)) and -Emm [fi] is the micromagnctic 
energy functional. In Eq. (JJJ r\ is a gaussian stochastic 
magnetic field with zero mean and correlations 



(? /CT (x, t)r]a> (x', 0) = <rS{t - t')a 3 S(x - x')S„ 



(2) 



where a is the (local) volume of the finite clement grid. 
The strength of the noise is given by a — 2ak-oT/h, pro- 
portional to the Gilbert damping parameter a and to 
the thermal energy k^T. (h is Planck's constant which 
relates energy and frequency.) In using this expression 
for the strength of the fluctuations we neglected the in- 
fluence of current on thermal magnetization fluctuations 
|2flj which is higher order j^l than the spin-torque effects 
studied here. 

The effective field can be separated into magnetic en- 
ergy and spin transfer torque contributions H = H| + 
H|j. For a ferromagnet with an easy z axis and a hard 
y axis we have 



hl = jv 2 n + 2cji« z - 2uj n v + h, 



(3) 



where J is the spin stiffness, H oxt is an external field 
in the easy-axis direction and u){ and lu q are respectively 
the easy axis and hard-plane anisotropy constants. As- 
suming only locality implied by smooth magnetization 
textures, the spin-transfer torque can be separated quite 
generally into contributions parallel and perpendicular to 
the spatial derivative of the magnetization: 



hi, x n = v. 



dr 



/3v s n x 



an 

Or 



(4) 



where the gradient is taken in the direction of current 
flow. In the absence of spin-orbit coupling, it follows 
from total spin conservation that (3 = and that the 
spin velocity v s = (J-\ — j±)/(—e(n-\ — nj.)), where j a and 
n a are majority and minority spin contributions to the 



2 



currents and spin-densities in the collinear limit. For 
realistic ferromagnets spin and orbital degrees of freedom 
are coupled, and the microscopic theory of i>s_and (3 is 
more challenging and still controversial [^laliaLillS]- The 
term proportional to (3, the nonadiabatic spin transfer- 
torque, plays a central role in the theory of current-driven 
domain wall motion. 

In the absence of current and noise, Eq. Q ad- 
mits time-independent solutions corresponding to do- 
main walls. In terms of the angles 9 and <j) defined 
by = (sin # cos 0, sin # sin 0, cosfl), the solution corre- 
sponding to an isolated domain wall centered at r& w is 
= and cos9(x — ra w ) = tanh[(x — rH W )/A], where 
the domain wall width A = yj ~j] (2wi) 3]. Rigid do- 
main wall motion is described by elevating rd w — > ^dw(^) 
and <fi(x,t) — > 4>o(t) to the role of collective dynami- 
cal variables 3J. The Langevin equations which describe 
their stochastic dynamics are 
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(5) 
(6) 



These equations can be derived heuristically by enforc- 
ing consistency with the stochastic LLG equations at the 
center of the domain wall. Alternately these equations 
can be derived by noting that the probability distribution 
P[Cl,t], generated by Eqs. and (J2J can be written as 
a path integral Pjtijl = / d[Q]5[£l ■ Cl - l] e - s ^, with 
effective action |3(X l3lL 132 



s[Ci] = 




x n-H 



(7) 



Inserting the domain-wall solution with time-dependent 
raw a- n d <po into this action gives rigid domain wall prob- 
abilities specified by the effective action 



.S'[,",u v . < = / dt' — 



t>o + a 



( - a4>o - w sin(20 o ) - y 



* 2 • 4 j 

-Wq sin 0o 
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(8) 



where N = 2AA/a 3 is the number of spins in a domain 
wall with cross-sectional area A. If we ignore the last 



term in this effective action, which vanishes in any case 
in the gaussian fluctuation limit, Eqs. (0 and © are 
recovered. This minor inconsistency in the theoretical 
treatment can be traced to the constant azimuthal angle 
across the domain wall in the variational ansatz. The 
final term must also be dropped if the rigid domain wall 
approximation action is to reproduce the correct equilib- 
rium probability distribution function. (Note that even 
when (3 = the Boltzmann equilibrium distribution is 
approached in the steady state.) The advantage of this 
approach is that we can read off the strengths of the 
gaussian noise terms r]^ and r\ r which are given by 



MW)) = (Vr(t)Vr(t')) = ~S(t - t') 



(9) 



Eqs. @i an d ® generalize the variational approach 
to current-induced domain- wall motion of Ref. 3] to fi- 
nite temperature. Tatara and Kohno Q derive their 
equations from an energy functional, an approach that 
has been criticized recently by Barnes and Maekawa |6|]. 
Our derivation does not appeal to an energy functional. 
Therefore, in addition to deriving the correct form of the 
noise terms required for a consistent treatment of thermal 
fluctuations, it also provides an alternative justification 
of the zero-temperature approach of Tatara and Kohno. 

These equations can be explored numerically without 
difficulty. However, it turns out to be possible to make 
some analytic progress. Solving for <po alone specializing 
to purely current-driven motion (H cxt = 0) we find that 



(l+a 2 )4> = -aw o sin(20o)H r -+r] < t > -ai] r . (10) 



This is the equation of motion for an overdamped Brow- 
nian particle which has been studied extensively with 
a large number of different physical motivations. Our 
theory for the equation of motion of the magnetization 
tilt should be contrasted with the treatment in Ref. [2l| 
which ultimately arrives at an underdamped Brownian 
particle equation to describe nonzero temperature do- 
main wall motion. 

The average (0o) can be calculated exactly for over- 
damped Brownian motion in a tilted periodic potential 
|33| . Inserting this well known result in Eq. (jSJ we find 
that 



0v B 27rAfc B T[l-e- 



^2-nHN(a-0)v s /{a\k u T) 



a hN{ e V(4>)/(k B T) d ^J^ e -VW)/(k B T) d ^^^^ e 2^HN( a -P)vJ( a Xk B T^ e -V (</,)/ {k B T) d(j) J* e V (<f,> ) / (k B T) d(j) , J 

(11) 
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FIG. 1: Domain wall velocity as a function of current, for 
various temperatures and values of j3. We take a = 0.02. The 
horizontal axis is normalized to v sc = \lu . 



where 

V(fa) = -hNcj o cos(2(b )/2 + HN(a - (3)v s cj) /(aX) , 

(12) 

is the tilted-washboard potential experienced by the 



Brownian "0o-particle" :3j. In the above equations we 
have assumed that a 2 <C 1, for notational convenience. 

It is well-known [3j|] that at zero temperature the equa- 
tion for an overdamped particle in the tilted periodic 
potential of Eq. I|10|) has solutions with (0) ^ only 
if | (a — (3)v s \ > a\uj . For (3 = the zero tempera- 
ture result for the average velocity of the domain wall is 
(^dw) oc y/ (v s /v sc ) 2 — 1 [3ll33j]. where the critical current 
for depinning the domain wall is v sc — Xu> 3] . For (3^0, 
the average domain wall velocity is nonzero at any finite 
value of v s even at zero temperature. 

In Fig. ^ we plot the average domain- wall velocity cal- 
culated from the full expression [Eq. (|llf> ] for various tem- 
peratures and values of (3. These curves were calculated 
for Gilbert damping parameter a — 0.02. We consider- 
only (3 < a since for (3 > a the domain wall velocity 
would decrease with increasing temperature in clear dis- 
agreement with experiment |l5| . One of the main conclu- 
sions of Ref. is that for (3 = and small currents the 
average domain-wall velocity (raw) °c exp(Cu s ) where C 
is a constant. In the low-temperature limit we find that 
for v s /v sc < 1 the average domain wall velocity is given 
by 



(rdw) = (Aw ) 2 - (v s f exp 

a v 




T-[— + — sin-M — 



sinh 



nhN (a - (3) v s 



(13) 



This difference in estimated domain wall velocities is 
closely analo gou s to the difference between the Langer- 
Ambegaokar |34l | and Halperin-McCumber estimates 
of phase slip resistance in thin superconducting wires. It 
follows from Eq. (|13fl that the observation of linear de- 
pendence of domain wall velocity on current does not 
necessarily imply that (3 ^ 0. A careful analysis of the 
temperature dependence of (rdw) will be necessary to de- 
termine the value of j3 from low current experiments, es- 
pecially so if Nhuj is not very large compared to ksT. 
Because of thermal noise, even when (3 = the domain 
is not intrinsically pinned at any nonzero temperature. 

Comparison with experiment — The results presented in 
Fig. n look qualitatively similar to the experimental re- 
sults of Yamanouchi et al. |l5| , in particular to the inset 
in Fig. 3 of Ref. [l^. A direct comparison with these 
experimental results is complicated by the fact that they 
are performed close to the critical temperature, with the 
result that the magnetic anisotropy energy-density and 
the polarization of the current depend on temperature. 
Moreover, the Gilbert-damping parameter a may also de- 
pend on temperature. The fact that curves at different 



temperatures in Fig. ^ are not as strongly offset verti- 
cally from each other at v s /v sc > 1 as the experimental 
results is most likely due to these additional temperature- 
dependent effects. 

One of the main results of Yamanouchi et al. is the 
empirical finding that for small currents ln(fd w ) cx ^/v^. 
From this the authors conclude that the domain wall un- 
dergoes current-induced creep motion at small currents. 
In Fig. |21 we plot the average domain wall velocity as a 
function of y/vl for various temperatures and values of 
(3. From this we observe that ln(fd w ) is proportional to 
the square root of the current over several decades of do- 
main wall velocity. It turns out that this approximate 
relationship is valid for a larger range of parameter space 
with decreasing temperature. In the experiments of Ya- 
manouchi et al. a typical velocity is (rdw) ~ 1 m s_1 - 
Putting this equal to Xuj we have for A = 17 nm and 
T ~ 100 K that k B T/{Nhw ) - 0.01 for N = 3.4 x 10 6 
Mn moments (we take the density of moments ./VMn ~ 1 
nm~ 3 ) in a wall of cross-sectional area A = 2 nm x 5 /im. 
For this value of k-B,T/(NhL0 o ) we conclude from Fig. [21 
that ln(fdw) oc y/v^ over a large regime, and that our 
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FIG. 2: Test of the scaling ln(rdw) oc 
peratures and values of j3. We take a - 



v s , for various tem- 
0.02. 



theory is therefore in agreement with the experimental 
results of Yamanouchi et al. 15]. Moreover, from Fig. [3 
we also conclude that for (3 — 0.01 this scaling does not 
hold. Hence, an additional conclusion from the experi- 
mental data is that j3 is much smaller than a in these 
ferromagnetic semiconductors. 

A similar estimate for typical parameters of the metal- 
lic nanowires used in studies of current-driven domain 
wall motion 0, 0, If3. Il3| leads to the conclusion that 
U-qT / (Nhw ) ~ 0.0002 and therefore that temperature 
effects on rigid domain wall motion are likely much less 
important in these systems. This difference arises mainly 
because the density of moments is roughly 40 times 
higher in metallic systems. 

In conclusion, we have presented a theory of the influ- 
ence of nonzero temperatures on current-driven motion 
of a rigid domain wall, and found qualitative agreement 
with ferromagnetic semiconductor experiments. An es- 
timate of the bending energy of the domain wall shows 
that these degrees of freedom are also thermally accessi- 
ble. It is, however, difficult to assess how they influence 
translation of the domain wall. The qualitative agree- 
ment of our results with the experiments of Yamanouchi 
et al. [l5j | indicates that thermal activation of rigid do- 
main wall motion could play an important role in these 
experiments. We expect that accounting for thermal fluc- 
tuations will also be important in assessing the impact of 
intended and unintended extrinsic domain wall pinning. 
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